Consequence of superfluidity on the expansion of a rotating Bose-Einstein condensate 
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We study the time evolution of a rotating condensate, that 
expands after being suddenly released from the confining trap, 
by solving the hydrodynamic equations of irrotational super- 
fluids. For slow initial rotation speeds, Qo, we find that the 
condensate's angular velocity increases rapidly to a maximum 
value and this is accompanied by a minimum in the deforma- 
tion of the condensate in the rotating plane. During the ex- 
pansion the sample makes a global rotation of approximately 
7r/2, where the exact value depends on Qq- This minimum 
deformation can serve as an easily detectable signature of su- 
perfiuidity in a Bose— Einstein condensate. 

PACS Numbers: 3.75.Fi, 67.40.Db, 67.90. +Z 



Recent experiments on trapped atomic Bose gases at 
very low temperature have revealed the existence of re- 
markable superfluid effects. Initial evidence of superfluid- 
ity came from the observation of quantized vortices [Qj^] 
and of the reduction of dissipative phenomena . 

Important signatures of superfluidity are also provided 
by the study of rotating condensates at low angular ve- 
locities, in the absence of quantized vortices, where the 
irrotationality of the flow implies the condition V x v = 
everywhere. In particular the quenching of the moment 
of inertia provides direct evidence of superfluidity, being 
related to the suppression of the transverse response of 
the many-body system Q . The recent observation of the 
scissors mode |^ has provided a clear example of this 
effect, confirming with high accuracy the predictions of 
theory §. The purpose of the present work is to investi- 
gate the consequences of superfiuidity on the expansion 
of an initially rotating and vortex-free condensate. 

In the absence of rotation the expansion of a con- 
densate after releasing the confining trap is well under- 
stood In particular the shape of a condensate, 
initially asymmetric in the x — y plane, when released 
approaches an asymptotic configuration where the role 
of the long and short axes is inverted with respect to 
the initial geometry. This inversion occurs because of 
the larger pressure felt by the system in the short direc- 
tion where the gradient of the density is higher. As a 
consequence the cross section of the condensate in the 
X — y plane will be instantaneously circular at some in- 
termediate time during the expansion. The time needed 
to reach this symmetrical configuration depends both on 
the initial deformation and on the presence of interatomic 
forces. For a noninteracting gas this time is given by 



{uxLOy)'^^'^ where tOx, iOy are the trapping frequencies of 
the harmonic potential which initially confines the gas. 
Interactions modify the time scale of the expansion in a 
profound way. In the Thomas-Fermi regime and in the 
limit of a very elongated condensate {u)y -C w^) the time 
needed to reach symmetry is given by ^l^y In this case 
the long axis does not expand significantly and the sym- 
metry of the configuration is achieved because of the fast 
expansion of the short axis j^]. 

If the condensate is initially rotating new features 
emerge during the expansion which are the object of the 
present investigation. In particular we predict that the 
condensate will never reach the symmetric configuration, 
because of the occurrence of a repulsive barrier caused 
by the irrotationality of the flow. 

An estimate of the effect is easily obtained by using 
angular momentum and energy conservation. Angular 
momentum conservation implies that when the system 
approaches the symmetric configuration in the x—y plane 
its angular velocity in the same plane becomes larger and 
larger as a consequence of the vanishing of the moment 
of inertia. This is a crucial consequence of the irrotation- 
ality of the superfluid motion. On the other hand the 
angular velocity cannot increase too much because of en- 
ergy conservation. By assuming that the initial energy of 
the system i?o is mostly converted into rotational energy, 
one finds Eq ^ fl^j.Qcr/'^- On the other hand angular 
momentum conservation requires f2o0o = ^cr^cr where 
ilo and Qq are the angular velocity and the moment of 
inertia at i = 0, while f2cr and Qcr are the corresponding 
values at the point of maximum rotational energy (here- 
after called critical time). One finds: Qcr ^ ^oQa/Eo 
and Qcr Eq/{QqQq). Since in a superfiuid the moment 
of inertia is given by the irrotational value Q — S^Qrig, 
where 
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is the deformation parameter of the cloud calculated 
with respect to its symmetry axes, and Orig — Nm < 
x^ -\- y^ > is the rigid value of the moment of inertia, 
one concludes that the system will reach a minimum de- 
formation which depends linearly on the initial angular 
velocity Qq, while its angular velocity, at the same time, 
will behave like I/Qq. 

The previous discussion suggests the following sce- 
nario: the rotating system will first expand in the short 
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direction (we consider here a condensate significantly de- 
formed in the x — y plane at t = 0). At a time, when the 
deformation nears its minimum value, the condensate will 
rotate rapidly causing the inversion of the long and short 
axes. The condensate then continues to expand along 
the original long axis. The occurrence of a minimum 
deformation during the exansion and the corresponding 
increase of the angular velocity can be regarded as a sig- 
nature of the superfluid behavior of the system. In fact 
a classical rotating gas expands in a very different way, 
because of the absence of the irrotationality constraint. 
In the classical case, neglecting the effects of collisions 
during the expansion, the angular velocity of the sample 
will decrease smoothly as a function of time according to 
the law riciit) = f7o/(l + f^gt^). 

In the following we will provide a quantitative descrip- 
tion of the expansion of a rotating Bose- Einstein conden- 
sate at zero temperature in the Thomas-Fermi regime. 
This regime is relevant to most current experimental situ- 
ations. The opposite limit, corresponding to the exapan- 
sion of a rotating ideal gas, will be considered in a sep- 
arate paper. In the Thomas-Fermi regime the equations 
of motion take the simplified form of the hydrodynamic 
equations of superfluids ||] . For the free expansion these 
equations, valid in the laboratory frame, can be written 
as an equation of continuity 
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V • (nv) = 0, 



and a Newton's 2nd-law "force" equation 
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where n and v are the condensate density and velocity 
field, respectively, and m is the atomic mass. The fac- 
tor g = ATrh^asc/m in the second equation measures the 
strength of the atom-atom interaction where Usc is the 
s-wave scattering length. 

The condensate is assumed to be held at equilibrium 
in a rotating trap just before turning off the external 
potential. In |10| it was shown that in the rotating frame 
the stationary solution of the equations of motion still 
corresponds to a spheroidal form for the density profile. 
At t = we can write 



and 



n (r, 0) = 6o - a-xox"^ - ciyoy'^ - dzoz^ 



V (r, 0) = aoV {xy) , 



(4) 
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where the x — y plane is perpendicular to the axis of 
rotation. The coefficients 6o, a.io and ao are determined 
by the number of atoms, TV, the trap frequencies, w^, 
and the initial angular velocity, fJo, and should fulfil the 
stationary criteria discussed in pc[ | . For example one has 
the relationship 



ao = 
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If the angular velocity fio is small, 6o and Oio coin- 
cide with their values in the absence of rotation andare 
given by 6o = and a,;o = "m^^fl'^g , where ^.tf is 

the Thomas-Fermi value of the chemical potential. In 
the same limit of small angular velocity the chemical 
potential fixes the initial radii, = (^o/oio)^^^ (see 
Eq. (^) of the condensate according to the relationships 
2^TF = JTiw^Xg = mUyYQ = uiloIZq. 

The solution of the expanding condensate can be writ- 
ten in the form 



n (r, t) =b (t) - ax (t) x - ay{t)y - {t) z 



- ay (t) y'^ - {t) 
a (t) xy (7) 



and 



V (r, t) = ^"^{ax (t) x^ + ay (t) y^ + a, (t) 

+ 2a (t) xy) . (8) 



The condensate occupies the region in space where n > 0. 
Note that the function b{t) defines the density at the 
center of the trap, n{0,t) ^ b{t). At t = one has 
6 (0) = bo, tti (0) = Qio, a (0) — ao and a (0) — ai (0) = 0. 

Equation (|^) shows that during the expansion the den- 
sity keeps a spheroidal form. Because of the initial an- 
gular velocity in the x ~ y plane the symmetry axis will 
rotate and the density will be no longer diagonal in the 
X, y coordinates. The angle of rotation can be easily 
evaluated in terms of the coefhcients of (|^) by intro- 
ducing the usual transformation x = x' cos (p + y' sin cj), 
y — —x'sm(j) + y'cos(f) and z = z' . This transforma- 
tion diagonalizcs the quantity = a^, (t) x^ + fly (i) y^ + 
ttz (t) z'^ + a {t) xy which takes the form 
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(a^ + ay){{x'y + {y'Y) + -{{a^ - ay) cos 20 



2' i^/vv / . V. / / . 2^ 

- asixi2<f){{x'f - {y'f) + a,{z'f 

with the angle of rotation fixed by the relation 

„ , a 
tan 20 = . 

^x 



(9) 



(10) 



The deformation parameter (|^), calculated in the rotat- 
ing frame, then takes the form (we assume 5 >Q) 



y/{ax - ay)^ + I 
(flx + a-y) 



(11) 



The equations of motion for the nine parameters a^, 
ai, a, a, and b are obtained by inserting Eqs. (0-^) into 
the hydrodynamic Eqs. (^) and (^. The equation of 
continuity yields 
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while the equation for the force provides the additional 
five equations 



Orr = 



— a„ = 
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2d + 2 (a^ + ay ) a — — ] a — 
d, + at - \ ^ ] a, =0 



(13) 



Equations (|lj) and complete the nine first-order or- 
dinary differential equations that we must solve in order 
to get n and v as a function of time. An equivalent for- 
malism was developed in [ pl| to investigate the behaviour 
of nonlinear oscillations carrying angular momentum. 

Equations (jlj) and (|3|) ensure the conservation of the 
integrals of motion, i.e., the number of atoms, N, the 
total energy E and the angular momentum, L^- These 
quantities can be expressed in terms of our variables. For 
example, the integral of the density over all space must 
equal N: 



N = 



Stt 
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while the angular momentum, given by the irrotational 

the form 

j2 + ay) 



law Lz — NmS^Qrig, takes the form 
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7 CLxO-y — a2/4 



(15) 



The conservation laws ( [14[ - |15| ) have been used to check 
the consistency of our numerical solutions. 
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FIG. 1. Angle of rotation of the condensate during the ex- 
pansion as a function of time for different values of the initial 
angular velocity. The choice of the parameters is given in the 
text. Time is given in dimensionless units (see text). 

It is important to note that the equations of motion 
(Eqs. (|lj) and (|l^)) are invariant with respect to the 
scaling transformation t — > t/A, ai — *■ Aa^, a — > Aa, 
tti — > A^ai,a — > A^a, so that the initial size of the con- 
densate can be actually excluded from the solution. For 
example, the dimensionless quantity S can be presented 
as S(tj, where t = t{2gaxQ/m)^/'^. For small initial an- 
gular velocities, t coincides with Lu^t where ujx is the fre- 
quency of the trapping potential. The function 6 (t) de- 
pends only on the initial shape of the condensate, i.e., on 
the ratios ayo/axo, azo/axo and on the dimensionless con- 
stant do = ao{2gaxo/m)~^^^ . The same behavior holds 
for the angle of rotation (f> (J) and for the dimensionless 

angular velocity Q (t) — d(f>/dt = Q{2gaxo/m)^^^^ . 

The results of the numerical integration of Eqs. ( |T^ ) 
and ([1^ ) are presented in Figs. |^|[ We have assumed 
that initially the condensate has a cigarlike shape elon- 
gated along the y-direction with azo/a-xo = ^,o-yo/axo = 
and A = 0.39. This corresponds to an initial deforma- 
tion in the x — y plane given by Sq — {1 — A^)/(l + A^) = 
0.74. Time is measured in units of the dimensionless 
quantity t. Three different values of the parameter Sq 
have been considered (0.05,0.15,0.45), corresponding to 
different choices of the initial angular velocity and of the 
trapping parameters. In the presence of rotation the re- 
lationships between the deformation parameters of the 
trap and the ones of the condensate are not trivial. 

Using the stationary solutions derived in [ p^ we find 
that the above choices for ao and A correspond to the 
values r^o = 0.07^2;, ujy — 0.39lUxi ~ ^-^^x'^ f^o — 
0.21072;, LOy — 0.451^2;, 10 z — ^ X ftud f^Q = ^ -2^0 x ^ oj y = 
1.41072, — 2.Qujx respectively. Notice that, except for 
values of SIq much smaller than the trapping frequencies 
LOxTiOy, the relationship between ao and f2o is not linear 
since the parameters axQ,ayQ entering Eq. (^) are renor- 
malized by the rotation of the trap. In particular the 
last choice = 1.207^, 07j, = 1.4107^) corresponds to 
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a configuration belonging to the overcritical branch dis- 
cussed in where the deformation of the condensate 
has opposite sign with respect to the one of the confining 
trap. 



condensate. 
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FIG. 2. Angular velocity of the condensate during the ex- 
pansion as a function of time for different values of the initial 
angular velocity (see Fig. 1.) 

Figure |l| shows that the angle of rotation (Tj reaches, 
for large times, a constant value. This is due to the 
rapid increase in the condensate size during the expan- 
sion which, in turn, causes a fast decrease of the angu- 
lar velocity. The actual value of the asymptotic angle 
depends on the initial conditions, differently from what 
happens in the ideal gas where (j) always approaches the 
value 7r/2. In Fig. || we show the angular velocity SI (Tj. 
One can see that, for small values of the initial angular 
velocity, fl exhibits a clear enhancement in the first stage 
of the expansion. Correspondingly the deformation pa- 
rameter S (Tj reaches a minimum (see Fig. ^ . If the ini- 
tial angular velocity is large, the minimum is shallow and 
eventually disappears. In this case also the angular ve- 
locity changes its behavior and decreases monotonically 
with time (see dotted line in Fig. H). In the recent exper- 
iment of Ref . , steady rotating states of magnetically 
trapped condensates were generated using rotating laser 
beams. The condensates were then imaged after expan- 
sion. By using the formalism of ref. the authors of 
Ref. [|l2| have checked that the deformation of the ex- 
panding condensate exhibits only minor changes. This 
result is the consequence of the relatively large values of 
the initial angular velocity considered in this experiment. 

In conclusion we have shown that the expansion of a 
rotating condensate reveals dramatic consequences of su- 
perfluidity associated with the reduction of the moment 
of inertia. These effects appear as a sudden increase of 
the angular velocity accompanied by the occurrence of a 
minimum in the deformation parameter. These effects, 
which are particularly pronounced for small initial an- 
gular velocities, should be easily observable by taking 
consecutive images during the expansion of the rotating 
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FIG. 3. Deformation parameter of an expanding conden- 
sate as a function of time for different values of the initial 
angular velocity (see Fig. 1). 

Useful discussions with Francesca Zambelli and David 
Feder are acknowledged. M.E. and C.W.C. acknowl- 
edge partial support from NSF grants nos. 9802547 and 
9803377. This work was supported by the Ministero della 
Ricerca Scientifica e Tecnologica (MURST). 



[1] M.R. Matthews, B.P. Anderson, P.C. Haljan, D.S. Hall, 
C.E. Wieman, and E.A. CorneU, Phys. Rev. Lett. 83, 2498 
(1999). 

[2] K.W. Madison, F. Chevy, W. Wohlleben, and J. Dalibard, 

Phys. Rev. Lett. 84, 806 (2000). 
[3] C. Raman, M. Kohl, R. Onofrio, D.S. Durfee, C.E. Kuk- 

lewich, Z. Hadzibabic, and W. Ketterle, Phys. Rev. Lett. 

83, 2502 (1999). 
[4] G. Baym, in Mathematical Methods in Solid State and Su- 

perfluid Theory, edited by R.C. Clark and E.H. Derrick 

(Oliver and Boyd, Edinburgh, 1969); S. Stringari, Phys. 

Rev. Lett. 76, 1405 (1996); F. Zambelh and S. Stringari, 

Phys Rev. A 63, 033602 (2001). 
[5] O.M. Marago, S.A. Hopkins, J. Arlt, E. Hodby, G. Hechen- 

blaikner, and C.J. Foot, Phys. Rev. Lett. 84, 2056 (2000). 
[6] D. Guery-Odelin and S. Stringari, Phys. Rev. Lett. 83, 

4452 (1999). 

[7] Y. Castin, R. Dum, Phys. Rev. Lett. 77, 5315 (1996). 

[8] F. Dalfovo, C. Minniti, S. Stringari, and L. Pitaevskii, 
Phys. Lett. A 227, 259 (1997). 

[9] S. Stringari, Phys. Rev. Lett. 77, 2360 (1996); F. Dal- 
fovo, S. Giorgini, L. Pitaevskii, and S. Stringari, Rev. Mod. 
Phys. 71, 463 (1999). 
[10] A. Recati, F. Zambelli, and S. Stringari, Phys. Rev. Lett. 
86, 377 (2001). 

[11] P. Storey and M. Ols hani, Phys. Rev. A 6 2, 033604 (2000). 
[12] K.W. Madison et al., |cond-mat/0101051. 



4 



